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We report an analysis of charmless hadronic decays of charged B mesons to the final state K'^-k^-k^ , 
using a data sample of 470.9 ± 2.8 million BB events collected with the BABAR detector at the T{AS) 
resonance. We observe an excess of signal events with a significance above 10 standard deviations 
including systematic uncertainties and measure the branching fraction to be B (B^ — >■ _R'"'"7r''7r'') = 
(15.5 ± 1.1 ± 1.6) X 10~®, where the uncertainties are statistical and systematic, respectively. 



Recent measurements of rates and asymmetries in 
B — > Kn decays have have generated considerable in- 
terest because of possible hints of new physics contribu- 
tions [ij, I4I. Unfortunately, hadronic uncertainties pre- 
vent a clear interpretation of these resuhs in terms of 
physics beyond the Standard Model (SM). A data driven 
approach, involving measurements of all observables in 
the B — )• Kt: system can in principle resolve the theoret- 
ical situation, but much more precise measurements (i.e. 
much larger data samples) w\\\ be needed [3|-|5|- 

It is interesting to study the related decays to 



pseudoscalar-vector final states B — > K*!: and B — > 
Kp p-|8|. In Table U we review the existing experimental 
measurements of the channels in the B — >■ K*^ system. It 
is evident that improved measurements of the K*+7r° [9| 
decay are needed. 

Due to the non-negligible width of the K* resonances, 
the quasi-two-body modes are best studied in the analy- 
sis of the three-body Dalitz plot. The four Ji'*7r decays 
populate six Kttii Dalitz plots (the four Kp decays also 
produce four of the same six final states) . Of these, Dalitz 
plot analyses oiK+TT+ir- [il,[l3l, KItt+tt- ^Wi and 



TABLE I: Experimental measurements of B 
Average values come from HFAG [id ]. 



K*TT decays. 



Mode 


Bx 10" 


Acp 


References 


K'^TT- 


10.3 ±1.1 


-0.23 ±0.08 


\11-U} 


K*+t:° 


6.9 ±2.3 


0.04 ±0.29 ±0.05 


[15] 


K*K+ 


go +0.8 
^■•^ -0.9 


-n 020 +°'"''' 

U.UZU _o.o61 


[16-i8] 


K*'\° 


2.4 ±0.7 


-0.15 ±0.12 ±0.02 


[12, i9] 



K~^TT~TT^ [12] have been performed to date. The first two 
of these have shown the presence of a poorly-understood 
structure, dubbed the /x(1300), in the tt+tt^ invariant 
mass distribution. A study of the K^tt'^tt'^ Dalitz plot 
would help to elucidate the nature of this peak, since 
even-spin states will populate both Ktt^tt^ and Ktt'^tt'^ 
(assuming isospin symmetry) , while odd-spin states can- 
not decay to Tr^Tr*^. 

Knowledge of the iiT+Tr^Tr'^ Dalitz plot may also help to 
clarify the interpretation of the inclusive time-dependent 
analyses ^ oi B° -)■ K^tt^tt^ ^^. Currently, these 
results show the largest deviation, albeit with a large 
uncertainty, among hadronic b ^>- s penguin-dominated 
decays [l^ from the naive Standard Model expectation 
that the time-dependent CP violation parameter should 
be given by Sep ~ —Vcp sin(2/3), where rjcp is the CP 

eigenvalue of the final state (+1 for K'^tt^tt'^) and B is an 

I — I nTi 

angle of the Cabibbo-Kobayashi-Maskawa [2J, [25| uni- 

tarity triangle. Such deviations could be caused by new 
physics, but in order to rule out the possibility of sizeable 
corrections to the Standard Model prediction, better un- 
derstanding of the population of the Ktt'^tt'^ Dalitz plots 
will be necessary. 

In this article, we present the results of a search for 
the three-body decay B^ — >■ K^tt^tt'^, including short- 
lived intermediate two-body modes that decay to this 
final state. This is the first step towards measuring the 
properties of contributing resonant modes. There is no 
existing previous measurement of the three-body branch- 
ing fraction, but several quasi-two-body modes that can 
decay to this final state have been seen, with varying sig- 
nificances. These include B~^ — » fo(980)K^, observed 
in the /o(980) — > tt+tt" channel [la, [l3l and also seen 
in /o(980) -^ K+IC^^, B+ ^ f2{l270)K+, seen in 
/2(1270)^7r^ 
in K*+{892) - 




and B+ -^ ii:*+ (892)71°, seen 
K+TT° (lii. The decay B+^^oK+ 
has also been observed with Xco 
K+K- 



M,^ 



M,m and 



XcO 

The data used in the analysis, collected with the BaBAR 
detector [28J at the PEP-II asymmetric energy e+e" col- 
lider at SLAC, consist of an integrated luminosity of 
429 fb~ recorded at the T(45') resonance ("on-peak") 
and 45fb~ collected 40 MeV below the resonance ("off- 
peak"). The on-peak data sample contains the whole 
BaBAR dataset of 470.9 ± 2.8 miUion BB events. 

We reconstruct i3+ — >■ X+7r°7r° decay candidates by 
combining a K^ candidate with two neutral pion candi- 



dates. The K~^ candidates are required to have a min- 
imum transverse momentum of 0.05 GeV/c and to be 
consistent with having originated from the interaction 
region. Separation of charged kaons from charged pi- 
ous is accomplished with energy-loss information from 
the tracking subdetectors, and the Cherenkov angle 
and number of photons measured by a ring-imaging 
Cherenkov detector. The efficiency for kaon selection 
is approximately 80 % including geometrical acceptance, 
while the probability of misidentification of pions as 
kaons is below 5 % up to a laboratory momentum of 
4 GeV/c. Neutral pion candidates are formed from pairs 
of photons with laboratory energies above 0.05 GeV and 
lateral moments between 0.01 and 0.6. We require that 
the mass of the reconstructed tt*' is within the range 
0.115 GeV/c^ < TO^^ < 0.150 GeV/c^ and that the ab- 
solute value of the cosine of the decay angle in the tt'^ 
rest frame is less than 0.9. We exclude candidates that 
are consistent with the B~^ — >■ K'^K'^, K'^ — )• tt^tt^ decay 
chain by rejecting events that contain a candidate that 
satisfies 0.40 GeV/ c^ < m^o^o < 0.55 GeV/ c^. This veto 
has a signal efficiency of at least 96 % for any charmless 
resonant decay and is almost 100 % efficient for nonreso- 
nant i?+ — > K+tt^tt^ and B+ — > XcoK^ decays. 

Due to the presence of two neutral pions in the final 
state, there is a significant probability for signal events to 
be misreconstructed, due to low momentum particles be- 
ing exchanged with particles from the decay of the other 
B meson in the event. We refer to these as "self-cross- 
feed" (SCF) events, as opposed to correctly reconstructed 
(CR) events. Using a classification based on Monte Carlo 
information, we find that in simulated events the SCF 
fraction depends strongly on the Dalitz plot distribution 
of the signal, and ranges from 2 % for i?+ -> XcoK~^ de- 
cays to 30% for B+ -^ f2{l270)K+ decays. 

To suppress the dominant background contribution, 
which arises from continuum e"'"e~ -^ qq {q — u, d, s, c) 
events, we employ a neural network that combines 
four variables commonly used to discriminate jet-like qq 
events from the more spherical BB events. These are 
the ratio of the second to the zeroth order momentum- 
weighted angular moment, the absolute value of the co- 
sine of the angle between the B direction and the beam 
axis, the absolute value of the cosine of the angle between 
the B thrust axis and the beam axis, and the absolute 
value of the output of a neural network used for "fiavour 
tagging" (i.e. for distinguishing B from B decays using 
inclusive properties of the decay of the other B meson 
in the T{AS) -^ BB decay ^). The first three quan- 
tities are calculated in the center-of-mass (CM) frame. 
We apply a loose criterion on the neural network output 
(NNout) which retains approximately 90% of the signal 
while rejecting approximately 82 % of the qq background. 

In addition to NNout , we distinguish signal from back- 
ground events using two kinematic variables: the dif- 
ference AE between the CM energy of the B candi- 



date and \fsl1-, and the beam-energy-substituted mass 
mES = \/s/4 — p^, where -y/s is the total CM energy and 
Pb is the momentum of the candidate B meson in the 
CM frame. The signal rriES distribution for CR events is 
approximately independent of the B'^ — ^ K'^-k^-k^ Dalitz 
plot distribution and peaks near the B mass with a res- 
olution of about 3MeV/c^. We select signal candidates 
that satisfy 5.260 GeV/c^ < mEs < 5.286 GeV/c^. The 
CR signal AS distribution peaks near zero, but has a 
resolution that depends on the signal Dalitz plot distri- 
bution, which is a 'priori unknown. To avoid possible bi- 
ases [30| we apply tighter selection criteria, — 0.15GeV < 
Ai5 < 0.05 GeV, which have an efficiency of about 
80% for signal while retaining only about 30% of the 
background (both compared to the looser requirement 
I A£:| < 0.30 GcV), and do not use AS in the fit described 
below. 

The efficiency for signal events to pass all the selection 
criteria is determined as a function of position in the 
Dalitz plot. Using a Monte Carlo (MC) simulation in 
which events uniformly populate the phase-space, we ob- 
tain an average efficiency of approximately 16 %, though 
values as low as 8 % are found near the corners of the 
Dalitz plot. 

An average of 1.3 B candidates is found per selected 
event. In events with multiple candidates we choose the 
one with the smallest value of a x^ variable formed from 
the sum of the x^ values of the two 7r° candidate masses. 
This procedure has been found to select the best recon- 
structed candidate more than 90 % of the time, and does 
not bias our fit variables. 

We study residual background contributions from BB 
events using MC simulations. It is found that these 
events can be combined into four categories based on 
their shapes in tties and AS. The first category com- 
prises two-body modes (mainly 5+ — > i^T+Tr"); the sec- 
ond contains three-body modes (mainly _B+ — ^ 7^*+7 
and B^ — ?> 7r+7r°7r°); the third and fourth are composed 
of higher multiplicity decays (many possible sources, 
with or without intermediate charmed states) with miss- 
ing particles, and are distinguished by the absence or 
presence of a peak in the tties distribution respectively. 
Based on the MC-derived efficiencies, total number of 
BB events, and known branching fractions [10|, |31|, we 
expect 70 ± 9, 39 ± 18, 1090 ± 40 and 170 ± 30 events in 
the four respective categories. 

To obtain the B^ — >■ iC+Tr^Tr^ signal yield, we perform 
an unbinned extended maximum likelihood fit to the can- 
didate events using two input variables: tties and NNout. 
For each component j (signal, qq background, and the 
four BB background categories) , we define a probability 
density function (PDF) 

P] = 7',(mEs^)7',(NNout^), (1) 

where i denotes the event index. The signal component 



is further separated into CR and SCF parts 

Vl,^ = (1 - /sCF)7'cR(mEs')7'cR(NNout') + 
/sCF^SCF(mEs07'sCF (NNout') 



(2) 



where /scf is the SCF fraction. The extended likelihood 
function is 



£=ne " n 



E-.^ 



(3) 



where n_,(fe) is the yield of the event category i{k). 

For the signal, the toes PDFs for CR and SCF are 
described by an asymmetric Gaussian with power-law 
tails and a third order Chebychev polynomial, respec- 
tively. Both CR and SCF NNout PDFs are described by 
one-dimensional histograms. We fix the shape param- 
eters to the values obtained from the B^ — > K^n^n^ 
phase-space MC sample, after adjusting them to account 
for possible differences between data and MC simula- 
tions determined with a high statistics control sample 



.0^ (^+^0^ 



decays. For the 



of B+ -^ D"p+ -^ {K^ 
continuum background, we use an ARGUS function 32 1 
to parameterize the toes shape. The continuum NNout 
shape is modelled with a parametric step function func- 
tion with 20 bins. One-dimensional histograms are used 
as nonparametric PDFs to represent all fit variables for 
the four BB background components. The free param- 
eters of our fit are the yields of signal and continuum 
background together with the parameters of the contin- 
uum Toes and NNout PDFs. All the yields and PDF 
shapes of the four BB background categories are fixed 
based on MC simulations. 

The results of the fit are highly sensitive to the value of 
/scf, which depends strongly on the Dalitz plot distribu- 
tion of signal events and cannot be determined directly 
from the fit. To circumvent this problem, we adopt an 
iterative procedure. We perform a fit with /scf fixed 
to an initial value. We then construct the signal Dalitz 
plot from the signal probabilities for each candidate event 
{sW eights) calculated with the sTlot technique 33[, and 
determine the corresponding average value of /scf- We 
then fit again with /scf fixed to the new value, and re- 
peat until the obtained values of the total signal yield 
(CR -f SCF) and /scf are unchanged between iterations. 
This method was validated using MC and was found to 
return values of /scf that are accurate to within 3 % of 
the nominal SCF fraction. Convergence is typically ob- 
tained within 3 iterations. 

We cross-check our analysis procedure using the high 
statistics control sample described above. We impose se- 
lection requirements on the D and p candidates' invari- 
ant masses: 1.84GeV/c^ < tox+tt-ttO < 1.88GeV/c^ and 
0.65GeV/c2 < to^+^o < 0.85GeV/c2. We fit the data 
with a likelihood function that includes components for 
the control channel, and for backgrounds from BB and 



qq backgrounds. We find a yield consistent within statis- 
tical uncertainties with the expectation based on world- 
average product branching fractions [3l[. 

We apply the fit method described above to the 31 673 
selected candidate B~^ —>■ ii'+Tr^Tr^ events. Convergence 
is obtained after four iterations with a yield of 1220 ± 85 
signal events and a SCF fraction of 9.7 %. The results of 
the fit are shown in Fig. [TJ The statistical significance of 
the signal yield, given by ^2 A In C where A In £ is the 
difference between the negative log likelihood obtained 
assuming zero signal events and that at its minimum, is 
15.6 standard deviations (a). Including systematic un- 
certainties (discussed below), the significance is above 
lOcr. 

The 5+ — > A'+7r''7r° branching fraction is determined 
from the result of the fit, dividing the signal sV^^ishts 
by event-by-event efficiencies that take the Dalitz-plot 
position dependence into account, and summing them to 
obtain an efficiency-corrected signal yield of 7427 ± 518 
events. We further correct for the effect of the K^ veto 
and a bias in the fitted signal yield, as determined from 
Monte Carlo pseudoexperiments generated with a sig- 
nal component with the same values of the yield and 
SCF fraction as found in the fit to data. Finally, we di- 
vide by the total number of BB events in the data sam- 
ple to obtain our measurement of the branching fraction 
B{B+ ^ K+Ti^TT^) = (15.5 ±1.1 ±1.6) X 10"^ where 
the first uncertainty is statistical and the second is sys- 
tematic. 

We assign systematic uncertainties due to (i) uncer- 
tainties in CR signal mEs PDF shapes (0.8 %) evaluated 
using the B+ -^ W p+ -^ {K+tt'-'k") (tt+t:°) control 
sample; (ii) uncertainties in CR signal and BB back- 
ground NNout PDF shapes (4.9 %) evaluated using uncer- 
tainties in the data/MC ratio determined from the con- 
trol sample and applying them to all PDFs in a correlated 
manner; (iii) uncertainties in the SCF signal mES and 
NNout PDF shapes (1.7% and 0.7%, respectively) evalu- 
ated considering a range of SCF shapes corresponding to 
different signal Dalitz plot distributions; (iv) uncertainty 
in the SCF fraction (2.5%) from varying the value used 
in the fit within a range of uncertainty determined from 
Monte Carlo pseudoexperiment tests of our iterative fit- 
ting procedure; (v) uncertainties in the BB background 
PDFs due to finite MC statistics (0.8%), determined by 
varying the contents of the bins of the histograms used to 
describe the PDFs within their errors; (vi) uncertainties 
in the BB background mES PDF shapes due to data/MC 
differences (1.6 %), evaluated by smearing the PDFs with 
a Gaussian with parameters determined from the control 
sample; (vii) uncertainties in the fixed BB background 
yields (1.4%), evaluated by varying these within their 
uncertainties; (viii) uncertainty in the correction due to 
fit bias (1.8%), which corresponds to half the correction 
combined in quadrature with its error; (ix) uncertainties 
in the efficiency, with contributions from tracking (0.4 %), 



TABLE II: Summary of systematic uncertainties. 



Source 


Uncertainty 


CR signal tties PDF 




0.8% 


CR signal and BB background NNout PDFs 


4.9% 


SCF signal tties PDF 




1.7% 


SCF signal NNout PDF 




0.7% 


SCF fraction 




2.5% 


BB background PDFs (MC 


statistics) 


0.8% 


BB background tties PDFs 




1.6% 


BB background yields 




1.4% 


Fit bias 




1.8% 


Subtotal 


6.5% 


Tracking efficiency 




0.4% 


Particle identification 




1.0% 


Neutral pion efficiency 




6.0% 


AE cut efficiency 




4.0% 


NNout cut efficiency 




3.0% 


K°s veto 




2.0% 


Nbb 




0.6% 


Total 


10.4 % 



kaon identification (1.0%), neutral pion reconstruction 
(3.0% per neutral pion, so 6.0% in total), AE (4.0%) 
and NNout (3.0%) selection requirements, the X!? veto 
correction (2.0%); (x) uncertainty in the number of BB 
pairs in the data sample (0.6%). The total systematic 
uncertainty on the branching fraction is 10.4%. Includ- 
ing only systematic uncertainties that affect the yield, 
the total is 6.5%. Table ITTl summarizes the systematic 
contributions. 

In summary, using the full BaBAR data sample of 
429 fb^^ collected at the T(4S') resonance, we observe 
charmless hadronic decays of charged B mesons to 
the final state K'^ir^ir'^. The signal has a signifi- 
cance above 10 cr, after taking systematic effects into 
account. We measure the branching fraction to be 
B {B+ -)■ K+Tr°7T°) = (15.5 ± 1.1 ± 1.6) x 10"^ This is 
the first step towards understanding the composition of 
the Dalitz plot of this decay and measuring the properties 
of contributing quasi- two-body modes. 
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FIG. 1: Projections of candidate events onto tties (left) and NNout (right) , following requirements on the other fit variable in 
order to enhance signal visibility. Points with error bars show the data, the solid (blue) curves represent the total fit result, 
the dashed (green) curves show the total background contribution, and the dotted (red) curves are the qq component. The 
dash-dotted curves represent the signal contribution. 
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